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In the Ginzburg-Landau approach, we describe a new phase in neutral two-flavor quark matter 
in which gluonic degrees of freedom play a crucial role. We call it a gluonic phase. In this phase 
gluonic dynamics cure a chromomagnetic instability in the 2SC solution and lead to spontaneous 
breakdown of the color gauge symmetry, the electromagnetic (7(1), and the rotational SO(3). In 
other words, the gluonic phase describes an anisotropic medium in which the color and electric 
superconductivities coexist. Because most of the initial symmetries in this system are spontaneously 
broken, its dynamics is very rich. 
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I. INTRODUCTION 

It is natural to expect that cold quark matter may exist in the interior of compact stars. This fact motivated 
intensive studies of this system over the past few years (for a review, see Ref. [l|). While these studies firmly 
established that the cold and dense quark matter is a color superconductor, they also revealed a remarkably rich 
phase structure in this system, consisting of many different phases. The question which phase is picked up by nature 
is still open. 

In this Letter we will describe a new phase in neutral and /3-equilibrated two-flavor quark matter. We call it a 
gluonic phase. The name reflects a crucial role of gluonic degrees of freedom in the structure of its ground state. More 
precisely, besides the usual color superconducting condensate of quarks, there exist (vector) condensates of gluons 
in this phase. These vector condensates cure a chromomagnetic instability in the two-flavor superconducting (2SC) 
solution |3| and lead to a dramatic rearrangement of the 2SC ground state, most notably, to spontaneous breakdown of 
SU{2)c X U{l)evi X 50(3)rot Symmetry down to S0{2),-ot- Here SU{2)c and [/(l)em are the color and electromagnetic 
gauge symmetries in the 2SC medium, and 50(3)rot is the rotational group (recall that in the 2SC solution the color 
5't/(3)c is broken down to SU{2)c). In other words, the gluonic phase describes an anisotropic medium in which 
the color and electric superconductivities coexist. As we will discuss below, there may exist a class of solutions with 
vector condensates of gluons. The solution described in this Letter is similar to that found in the gauged cr-model 
with a chemical potential for hypercharge although physics in the present system is much richer. In particular, as 
will be shown in Sec. exotic hadronic states play an important role in its dynamics. 

The framework of the present study is the Ginzburg-Landau (GL) approach. The basic point in our analysis is the 
inclusion of light gluonic degrees of freedom in the GL effective action. This leads us to revealing the gluonic phase 
in neutral two-flavor quark matter. 

II. EFFECTIVE POTENTIAL WITH VECTOR CONDENSATES 

We study dense two-flavor quark matter in /3-equilibrium. For our purpose, it is convenient to use a phenomenolog- 
ical Nambu-Jona-Lasinio (NJL) model, more precisely, a gauged NJL model including gluons. Although usually the 
NJL model is regarded as a low-energy effective theory in which massive gluons are integrated out, we introduce glu- 
onic degrees of freedom because the gluons of the unbroken SU{2)c subgroup of the color SU{i)c are left as massless. 
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and, under certain conditions considered below, some other gluons can be also very light. As discussed in Ref. 
the confinement scale Aq^q in the two-flavor color 2SC phase is estimated as 0(10 MeV) or smaller. Thus it is not 
peculiar to consider gluonic degrees of freedom at energies less than the quark gap A. 

We neglect the current quark masses and, as usual in studying the 2SC phase, assume that the color superconducting 
condensate does not break parity. The Lagrangian density is then given by 



where 



(g'^i£e°75g)(gi£e°75'7^) 



^9, (1) 



C, = -\f^,F^^'^ (2) 

and 

i^^ ^ - ^5A^T^ F^^^d^AZ-d^A'^^+gpf^^AlAl. (3) 

Here A'^ are gluon fields, T" are the SU{i) matrices in the fundamental representation, and and t"^"^ are the 
antisymmetric tensors in the flavor and color spaces, respectively. In /3-equilibrium, the elements of the diagonal 
chemical potential matrix jl for up (u) and down [d) quarks are 

Mur = = A* - Hdr = IJ-dg = jj- + (5^^, (4a) 

tiub = jj-- 5^L- ^is, fj.db = P' + Sn- Us, (4b) 

with 

M = M-Y + y> <5M=Y- (5) 

Here the subscripts r, g, and b correspond to red, green and blue quark colors, /i is the quark chemical potential (the 
baryon chemical potential is given by /is = 3/i), fie is the chemical potential for the electric charge, and fis is 
the color chemical potential. The latter is simply connected with the vacuum expectation value (VEV) of the time 
component of the 8-th gluon : 

giO^^m- (6) 

By using the auxiliary field A" ~ icp ee°'^^q, the Lagrangian density can be rewritten as 

£ = -qiilf) + /i7°)9 - iA"[zg-£e"75<Z^] - \^i<f e,'^ ^^qW^ _ ^ + (7) 



We now introduce the Nambu-Gor'kov spinor. 



* = ( „c ) ■ (8) 



The inverse propagator of ^ including gluons is written as 

^ "I A+ [G„-J-ij' 

with 

[G+ ^ (/ + A - '5/xr3 - ii^U)^^ - 7 • P + arT'^. (10) 
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[Gq- ^ (p" - /i + SfiT^ + Ai8l,)7" - 7 • p - g^T"^ (11) 

and 

A- = -ie''e-/5A, A+ = -ie^'e^/^A* . (12) 

Here the diquark condensate is chosen along the blue color direction, the constant fields A"^ represent possible vector 
gluonic condensates in the model, and we introduced matrices t'^ = diag(l,— 1) and 1^ = diag(0,0, 1) acting in the 
flavor and color spaces, respectively. The effective potential in this model includes both giuons and the scalar field 
A. It is: 

^''^ = ^ + ^r^^r'^K^z^''^^" -\j -^i^det^;!. (13) 

We will utilize the hard dense loop approximation, in which only the dominant one-loop quark contribution is taken 
into account, while the contribution of gluon loops is neglected. On the other hand, we keep the tree contribution of 
giuons in the effective potential (|13|l . This is because we want to compare this contribution with that of hard dense 
loops in order to check the consistency of the hard dense loop approximation. 



III. DYNAMICS OF GLUONS IN THE GINZBURG-LANDAU APPROACH 



Let us consider dynamics of light gluonic degrees of freedom in this model. The two point function of giuons can be 
calculated from Lagrangian density {Tj). In Refs. '5, "6], the Debye and Meissner screening masses of the giuons in the 
2SC phase were explicitly calculated. For the giuons of the unbroken S'C/(2)c, i.e., A'^^\ A^'^\ and A'~^\ both Debye 
and Meissner masses vanish in the region (5/i < A (henceforth, for clarity, we will put color indices of gluon fields in 
parentheses). For the giuons yl^'*)"'^^^ the Meissner mass is approximately 



1 - ^ ' '^M < A. (14) 



Thus, near the critical point (5/i = A/\/2, the Meissner mass for A'^^^ is very small. As (5/x exceeds the value A/\/2, 
becomes negative, thus signalizing a chromomagnetic instability in the 2SC solution 0. On the other hand. 
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"''M,4 

around the critical point 5/i = A/V2, the SU{2)c singlet gluon A^^'' is heavy. Actually, it has nonvanishing Debye 
and Meissner screening masses in the whole region (5/i < A. This fact allows us to pick up the giuons A^^'>~^'^'> as 
relevant light degrees of freedom in the low energy effective theory around the critical point (5/i = A/^/2. Our goal is 
to describe the dynamics near this critical point in the GL approach. 

The quark gap A 7^ breaks the QCD symmetry SU{3)c down to SU{2)c. With respect to SU{2)c, the adjoint 
representation of SU (3)c is decomposed as 

8 = 3©2®2®1, (15) 

i.e., 

{4")} = {aW , 42) ^(3) ) e e e 4«) , (a = i, 2, • • • , 8). (16) 

Here we defined the complex doublet of the "matter" field describing color vector "kaons" : 
We also define 

m ^ 9^4) - 9.4) + <?e'™"4™)4"), (l, m, n = 1, 2, 3) (18) 
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and 

V,^^,-^gA<^^^. (19) 

Then the building blocks of the effective action are the following six, 

A'o, K,, Vo, V„ /oj, /,fe. (20) 

The effective action should be of course invariant under all initial symmetries, in particular, under the color SU (2)c 
and the rotational SO{'i)mt- 

Because of the chromomagnetic instability at (5/i > A/\/2, it is natural to study a spontaneous breakdown of the 
SU{2)c via the formation of a vector condensate {Kfj_) ^ 0. Because the instability is chromomagnetic, we assume 
that a spatial component of has a VEV. By using the rotational symmetry S'0(3)rot, one can take (K^) ^ 0. And 

because of the SU{2)c symmetry, without loss of generality, we can choose {Af^) ^ 0. This VEV breaks both the 
SU{2)c and the SO{i),ot- 

The following remarks are in order. The complex doublet K-^ plays here the role of a Higgs field responsible 
for spontaneous breakdown of the SU{2)c. The situation is similar to that taking place in the electroweak theory. 
The essential difference however is that now the Higgs field is a spatial component of the vector field leading also 
to spontaneous breakdown of the rotational symmetry. In this paper, the unitary gauge will be used in which 
K'^ — :^(0, (4^'') +03^'') where the real field af^ describes quantum fluctuations. The important point is that in the 
unitary gauge all auxiliary (gauge dependent) degrees of freedom are removed. Therefore in this gauge the vacuum 
expectation values of vector fields are well-defined physical quantities. 

With a broken SU{2)c, the SU (2)c gluons could have VEVs. A similar situation takes place in the gauged cr-model 
with a chemical potential for hypercharge, where the gauge SU{2)l is broken Motivating by that model, we 
assume 

{4'^)^0, (21) 

and use the following notation, 

i3^.g(Af), C^5(4'^ ^^.9(4'^) (22) 

(note that g{Af^ ) can be considered as a chemical potential /is related to the third component of the color isospin). 
As will be shown in Sec. IIVI such a solution with nonzero B, C and D vector condensates exists in this model indeed. 

We now describe a general symmetry breaking structure in this model. In the presence of /ie, the chiral symmetry 
is explicitly broken. Then the initial symmetry is 

[^C/(3)c]local X [C/(l)e™ X [/(1),3 X C/ ( 1 ) ,3 ] global X 50(3),ot, (23) 

where J7(l)^3 ^ are J7(l)-part of the chiral symmetry SU{2)l ji (note that there is no photon field in the model). The 
baryon charge is incorporated in the subgroup, 

S= il/®lc = 2(Q-/3), (24) 

where Q = diag(2/3, — 1/3) and /a = diag(l/2, — 1/2) acting on the flavor space. The quark gap A breaks both S'/7(3)c 
(down to SU{2)c) and U{l)em but a linear combination of and Q remains unbroken. The new electric charge of 
the unbroken U{l)em is 

Q = g--^r«. (25) 

The baryon charge is also changed to 

^ = 2(g-/3). (26) 
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The VEV (Ag ') breaks SU(2)c, but a linear combination of the generator from the SU{2)c and Q, 

Q^Q-T^^Q- - (27) 

v3 



determines the unbroken U{l)em (the new baryon charge is B = 2{Q — However, because does not commute 

4^ 



with T'^, the VEV (^3^'') breaks f/em(l)- The baryon charge is also broken. 



After all, we have: 

[5f/(3)Jlocal X [(7(l)e™ X U {lU X (7(1), 3] global X SO{3),ot 



^ [5C/(2),]iocal X [Uil)em X U (lU X C/(l),3 ]global X 50(3)rot (28) 



[[/(l)e™ X [7(1)^3 X ;7(l),3]giobal X SO{2)rot (29) 
[f/(l),3 X [/(l),3]giobal X 50(2)rot. (30) 



Thus, this system describes an anisotropic medium in which both the color and electric superconductivities coexist. 

Let us apply the GL approach to this system near the critical point ~ A/-\/2. The SU{2)c and S'0(3)rot 
symmetries dictate that the general GL effective potential, made from building blocks H2U|) and including operators 
up to the mass dimension four, is 

Ves = Va + \mIB^ + TDB'' + IascS'C^' + \\bdB^D'' + ]^\cdC^D'' + \\bB^, (31) 

where Va is the 2SC part of the effective potential. Here, while the coefficients As, XbCi ^bd, and Xcd are 
dimensionless, the dimension (in mass units) of the coefficient T in the triple vertex is one. Expanding the potential 
(|13|l with respect to B, C, and D, we can determine the potential H31|l . Amazingly, the structure of the gluonic part of 
effective potential (|31(l is quite similar to that of the potential in the gauged tr-model with the hypercharge chemical 
potential although the present system is much richer. 

Before realizing explicit calculations, we clarify the behavior of the effective potential (|31|) near the critical point. 
The stationary point of the effective potential H31() is given by equations 



dVos 
dB 

dC 

dD 

and 



= B[mI + \bB'^ + 2TD + \bcC^ + XbdD^]=Q, (32) 
= C[\bcB^ + XcdD^] =0, (33) 
= TB^ + \bdDB^ + XcdC^D ^Q, (34) 



^''^ 

We can expand /^e, MS: and A around B = C = D = 0, 

He = Me + Ce, (36) 

fis = Ms + Cs, (37) 
A = A + a, (38) 

where the bar-quantities are the 2SC solution, when B = C — D — 0. 

Let us assume that the origin (bifurcation point) of the solution with nonzero B, C, and D corresponds to a 
second order phase transition (as will become clear in a moment, this assumption is self-consistent). Then, taking an 
infinitesimally small B near the critical point, we easily find that 

ee,6,eA~0(i?2). (39) 
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It implies that the difference of Va in the new solution and that in the 2SC one is 

VAiA^°\^^T\^^f) - Va{a, fi,,fis) - o{b^). (40) 

This fact will be useful in our analysis below. 

From Eqs. H32|) - (|34|l . we find that when the 2SC solution becomes unstable (M^ < 0), a new solution occurs, if 
the parameters \bc and \cd satisfy 

Xbc > 0, XcD < (41) 
(in the next section, it will be shown that this constraint is satisfied indeed). The new solution is: 



^^°> = WV^' ^^°'"V3A^' ^^^^ 

where we neglected higher order terms of Mg. In Eq. (|42|l the conventions B > and C > are chosen. 
Near the critical point = 0, the solution behaves as 



S,oioc-M|, Csoicx^-M|, Dsoioc-A/|. (43) 

These scaling relations are quite remarkable. While the scaling relation for C is of engineering type, those for B and 
D are not (the origin of this is of course in the presence of the dimensional coefficient T in Eq. H42|l ). Such a scaling 
behavior implies that the B** and B^D^ terms in the effective potential are irrelevant near the critical point Mg — 0. 
Omitting them, we arrive at the reduced effective potential: 

Ves ^Va + ^MlB^ + TDB'^ + Iabc^'C^ + ]^\cdC^D''. (44) 

Notice that (V^ff — Va) ^ 0{B^). This fact and Eq. H40() imply that in the leading approximation one can use the 
bar-quantities, defined in Eqs. 1)36(1 - H38|l . in calculating Va, Mg, T, Xbc, and Xcd in the reduced potential. In other 
words, the effective potential can be decomposed into the "constant" 2SC part Va, with frozen fermion parameters, 
and the dynamical gluonic part: 

V,s ^ V^ff (A, Ae, Ms; B, C, D) = Fa (A, ^e, Ms) + \mIB'' + TDB"" + ^XbcB^C^ + ^XcdC^D^. (45) 
Then Eq. (|42|l is the exact solution of the potential H45|) and the energy density at the stationary point is found as 

Feff (A, fie, Ms; Ssol, CsoU I^sol) = ^^A + ^MIBI, ^ Va " ( ) < ^^A- (46) 

Therefore the gluonic vacuum is more stable than the 2SC one. 
By using the Gauss's law constraint 

TB^ + XcdC^D^O, (47) 

(3) 

we find the true effective potential without the non-dynamical degree of freedom Aq : 

V;Gau.s ^y^^ 1^^252 ^ ^XbcB'C - (48) 

It is easy to show that solution H42|l is a minimum by analyzing the curvature of T/^'^""**. 

In the next section, we will calculate Mg, T, Xbc, and Xcd- In particular, it will be shown that constraint (|41|l is 
satisfied near the critical point. 
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IV. DYNAMICS IN ONE-LOOP APPROXIMATION 



In this section, we determine the GL effective potential (I45|) in one-loop approximation and determine the dispersion 
relations for quarks in the gluonic phase. The 2SC Va part of the potential is known 



yA(A,/ie,^8) 



Me 



4G'a 127r2 



j4b_ 
127r2 



37r2 



A2 
7r2 



-A' 



4A2 



In-r^ 



A2 



A2 

2 



TT 



;a2 



is^l < A). 



(49) 



Here A is the ultraviolet cutoff in the NJL model and fiub, l^db, and ji are given in Eqs. Q and (|SJ). For clarity of 
the presentation, the bars in A, /ig and /is were omitted {0{ji^ / h?') and C'(A2/A2) and higher terms are neglected in 
this expression). Note that the color and electrical charge neutrality conditions in the 2SC solution yield [3| 



and 



A 3 1 

10 5 



2A 



(50) 



(51) 



which is consistent with the result of Ref. [3, ~ 0{A? / in the case of 5^ — 0. The size of A is essentially 
determined by tuning the NJL coupling constant Ga and cutoff A. 

After straightforward but tedious calculations of relevant one-loop diagrams from the fermion determinant in 
Eq. \i'6\ . we find the following relations in the region < A: 



Ml 



9' 




A2 


1 


~4 


-1 + 




807r2 


A2 





A2 



"^CD — 7 

9^ 



T 



1 

187r2 A2 



Ms M^ 



Ms 

2^2 ^ 247r2 A2 



M 

A4 



487r2 



A4 



Here the tree contribution of gluons 



1 



2 oi oi 



(52) 
(53) 
(54) 
(55) 

(56) 



was also taken into account. 

We see that the coefficient Xcd is definitely negative. The parameter Mg, which is expressed through the Meissner 
mass (|T^. is negative when 



(5/i > SfXc 



A I 3tt III 



V2 



2as p 



El 

Air' 



(57) 



Relation ((SUJl and Eq. ^ yield 



^=10^ 



rMs, 



(58) 
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and, at the critical point, we find from Eqs. H50|l . (|51|l and (|57() . H58() that /is is approximately 



In 



200Af 



12 



In 



200A2 O 



For realistic values A = (1.5 — 2.0)fi and as = 0.75 — 1.0, we obtain numerically 

^4-0-03-0.1. 



(59) 



(60) 



This implies that the tree gluon contribution decreases the value of 6fici- by 1.5%-5% in comparison to its value in the 
(non-gauged) NJL model. The smallness of this correction is in accordance with the dominance of hard-dense-loop 
diagrams. 

Let us now turn to the coefficient Xbc H53|) . At the critical point Sfi = (5/icr, it is: 



807r2 A2 



(61) 



Because the /ig//i'*-term is negligibly small, we conclude that the coefficient Xbc is positive near the critical point. 
Thus, constraint (|41ll is satisfied indeed. 

Utilizing Eqs. I|52(l - (|55|1 in Eq. (|42|l . one can obtain the solutions for B, C, and D in the near-critical region. Indeed, 
neglecting higher order terms in /^s/m in l|52|l - (|55|l . we get the approximate relations 



67r2 



1 - 



Xbc 



1607r2 



A 



1 



CD 



1 



Ana, 47r2 ' 



T ■ 



Ms 



167r2 167r2 



-), (62) 



which lead us to the near-critical solution: 



(5/i2 — 5^ 



2 16^1/10 1 



27 



1 + f (3 



2ii 



(63) 



sol 



^JSfj? — (5/i2j. 4:^/5 p. 
(5/i2 - Sn'^^. I 



n 




(3H 


2ir \ 








as J . 



(64) 
(65) 



It is noticeable that this solution describes nonzero field strengths F^"' which correspond to the presence of non 



abelian constant chromoelectric-like condensates in the ground state: 



E. 



E. 



(2) 
3 

(7) 



(2) 



1 



03 



Cso\D 



F, 



(7) 



03 



2g 



-Bsol ( 2/i8 — Dso\ ) ■ 



(66) 
(67) 



We emphasize that while an abelian constant electric field in different media always leads to an instability, ^ non- 
abelian constant chromoelectric fields do not in many cases: For a thorough discussion of the stability problem for 



^ In metallic and superconducting media, such an instability is classical in its origin. In semiconductors and insulators, this instability is 
manifested in an creation of electron-hole pairs through a quantum tunneling process. 
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constant SU{2) non-abelian fields in theories with zero baryon density, see Ref. i. On a technical side, this difference 
is connected with that while a vector potential corresponding to a constant abelian electric field depends on spatial 
and/or time coordinates, a constant non-abelian chromoelectric field is expressed through constant vector potentials, 
as takes place in our case, and therefore momentum and energy are good quantum numbers in the latter. 

In order to illustrate the stability issue in the gluonic phase, let us consider the dispersion relations for quarks there. 
Because the vacuum expectation values H63(l -()65 |l are small near the critical point and because red and green quarks 
are gapped in the 2SC phase, the dispersion relations for gapless blue up and down quarks are of the most interest. 
From Eq. @ we find that up to the first order in they are 

The B^-terms in Eqs.(EH|) and (EU lead to non-spherical Fermi surfaces determined by the following equations: 

1-1 -S^isin^e 1 S^icos^e Me + MS X 

= Mu6 —. ■ ^2 -rj— — ^, (blue up 70 

P^^fMb-^^. + ^^A A2 ^ 2 ^ blue down 71 

where we neglected higher order terms of B^ and defined the angle 9, 

p^ = \p\ cos 9. (72) 

The dispersion relations H()8|) and (|69|l clearly show that there is no instability in the quark sector in this problem. 

As to bosonic degrees of freedom (giuons and composite bosons) , because it is very involved to derive their derivative 
terms from the fermion loop in the gluonic phase, this issue is beyond the scope of this letter. It is however noticeable 
that there are no instabilities for bosons in a phase with vector condensates in the gauged cr-model with a chemical 
potential for hypercharge |3|. Although that model is much simpler than the present one, its phase with vector 
condensates has many common features with the gluonic phase and this fact is encouraging. 

We emphasize that these constant color condensates in the gluonic phase do not produce long range color forces 
acting on quasiparticles. This can be seen from the dispersion relations (68) and (69) for quarks in this model. They 
show that momentum and energy are conserved numbers. It would be of course impossible in the presence of long 
range forces. The role of these condensates is actually more dramatic: They change the structure of the ground 
state, making it anisotropic and (electrically) superconducting. Only in this sense, one can speak about a long range 
character of the condensates. 



V. MORE ABOUT DYNAMICS IN THE GLUONIC PHASE 



In this section, we will describe some additional features of the gluonic phase. In particular, we will point out that 
a condensation of exotic vector mesons takes place in this dense medium. 

The gluonic vector condensates are mostly generated at energy scales between the confinement scale in the 2SC 
state, which is ^ 10 MeV, and the baryon chemical potential, which is about 300-500 MeV. It is the same region 
where the chromomagnetic instability in the 2SC phase is created and where the hard dense loop approximation is 
(at least qualitatively) reliable. At such scales, giuons are still appropriate dynamical degrees of freedom and utilizing 
the Higgs approach with color condensates in a particular gauge is appropriate and consistent: It is a region of hard 
physics. Because the gluonic phase occurs as a result of a conventional second order phase transition, the vector 
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condensates are very small only in the immediate surroundings of the critical point 6fi = A/^/2. Outside that region, 
their values should be of the order of the typical scale (5/x ~ A ~ 100 MeV. 

These condensates represent hard dynamics connected with the appearance of a new parameter, chemical potential 
for the electric charge fie ^ 100 MeV. As a result, the SU{2)c gauge symmetry becomes completely broken and the 
strong coupling (confinement) dynamics presented in the 2SC solution at the scale of order 10 MeV is washed out. 
In other words, a conventional Higgs mechanism is realized in the gluonic phase. In this respect, the gluonic phase 
is similar to the color-flavor locked (CFL) phase, where the constant color condensates (although not vector ones) 
completely break the SU{3)c color gauge symmetry 

It is easy to check that the electric charge Q and the baryon number B introduced in Section 3 are integer both 
for gluons and quarks. Do they describe hadronic-like excitations? We believe that the answer is "yes". The point 
is that in models like this one, with a Higgs field in the fundamental representation of the gauge group, there is no 
phase transition between Higgs and confinement phases |^. These two phases provide dual, and physically equivalent, 
descriptions of dynamics (the complementarity principle). In particular, they provide two complementary descriptions 
of a spontaneous breakdown of global symmetries, such as the rotational S0(3) and the electromagnetic U{1) in the 
present case. Following Ref. , one can apply the dual gauge invariant approach in this model and show that all the 
gluonic and quark fields can be replaced by colorless composite fields. The flavor numbers of these fields are described 

by the conventional electric and baryon charges Q and B. They are integer and coincide with those the operators Q 

and B yield for gluonic and quark fields. 

While these issues will be considered in more detail elsewhere, here we would like to point out the following 

noticeable feature of these states: some of them are exotic. For example, the electric and baryon charges Q and B of 
A'i^'' = A^^^ + zA^^^ gluons are equal to +1 and +2, respectively. Because Ag^^ gluons are condensed in the gluonic 
phase, we conclude that in the dual gauge invariant description this corresponds to a condensation of exotic vector 
mesons. In this regard, it is appropriate to mention that some authors speculated about a possibility of a condensation 
of vector p mesons in dense baryon matter 10] . The dynamics in the gluonic phase yield a scenario even with a more 
unexpected condensation. 

VI. SUMMARY AND DISCUSSIONS 

The gluonic phase whose existence was shown in this paper is very different from all known phases in dense quark 
matter discussed in the literature. Also, to the best of our knowledge, no phase like that has been considered in 
condensed matter physics. One of its features is the presence of non-abelian constant chromoelectric condensates in 
the ground state. They make the dynamics of the gluonic phase to be manifestly non-abelian. ^ 

Because most of the initial symmetries in this system (including the rotational S'O (3)rot and the electromagnetic 
U{1)) are spontaneously broken, the spectrum of excitations in the gluonic phase should be very rich. In particular, 
there should be two gapless Nambu-Goldstone (NG) modes connected with the two broken generators of the rotational 
group and one NG mode corresponding to the broken electric charge (the latter mode will be absorbed into photon field 
neglected in our model). Another interesting feature of the gluonic phase is that there are excitations corresponding 
to exotic hadrons. We are planning to return to this issue elsewhere. 

The solution Ht)3|l - Ht)5() described in this paper corresponds to a minimum of the effective potential. Whether or 
not this minimum is global is an open question. We suspect that there may exist a class of solutions with vector 
condensates. In this regard, it is instructive to describe the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase 
from this point of view. It is easy to show that the LOFF solution with one plane wave along, say, third spatial 
coordinate can be gauge transformed into a solution with a usual (homogeneous) diquark condensate and a vector 



^ The role of constant chromomagnetic non-abelian fields in the dynamics of color superconductivity was studied in Ref. Illl . It was shown 
that they could enhance the value of a diquark condensate. Unl ike the gluonic phase, where non-abelian chromoelectric condensates are 
solutions of dynamical equations, the chromomagnetic fields in Illl are external. 

^ The spectrum of excitations in the gauged cr-model with hypercharge chemical potential |3, where vector condensates also occur, strongly 
supports these expectations. 
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condensate (Aij '). [Note that there is no such a transformation in the case of the LOFF solution with two or more 
plane waves |13|.] Unlike the present gluonic solution, there are no non-abelian chromoelectric condensates in that 
case, and, therefore, the LOFF dynamics is not genuinely non-abelian. Still, because the LOFF solution has been 
recently used to cure a chromomagnetic instability in the neutral 2SC phase Il3|, it would be worth to study a 
possibility of extending the present solution by including the VEV {Af''^). 

Last but not least, it would be interesting to check the possibility of the existence of a gluonic phase in three- flavor 
dense quark matter. It is especially interesting because a chromomagnetic instability has been recently revealed also 
in that case [l6l |. 
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* A model for describing dynamics for larger ~ A, close to the edge of the gapless phase 0, has been recently considered in Ref. 
It is an abelian gauge model with a photon gauge field having a constant VEV. By using an appropriate gauge transformation, one can 
remove this VEV, with a cost of introducing a phase factor in the order parameter, as it takes place in the LOFF state. 



